Stability of a Pivoting Strategy for Parallel Gaussian Elimination
نویسنده
چکیده
Gaussian elimination with partial pivoting achieved by adding the pivot row to the kth row at step k, was introduced by Onaga and Takechi in 1986 as a means for reducing communications in parallel implementations. In this paper it is shown that the growth factor of this partial pivoting algorithm is bounded above by μn < 3 n−1, as compared to 2n−1 for the standard partial pivoting. This bound μn, close to 3 n−2, is attainable for a class of near-singular matrices. Moreover, for the same matrices the growth factor is small under partial pivoting. AMS subject classification: 65F05, 65G05.
منابع مشابه
Complete pivoting strategy for the $IUL$ preconditioner obtained from Backward Factored APproximate INVerse process
In this paper, we use a complete pivoting strategy to compute the IUL preconditioner obtained as the by-product of the Backward Factored APproximate INVerse process. This pivoting is based on the complete pivoting strategy of the Backward IJK version of Gaussian Elimination process. There is a parameter $alpha$ to control the complete pivoting process. We have studied the effect of dif...
متن کاملThreshold-pivoting in parallel Gaussian elimination for improved efficiency
The use of threshold pivoting with the purpose to reduce fill-in during sparse Gaussian elimination has been generally acknowledged. Here we describe the application of threshold pivoting in dense Gaussian elimination for improving the performance of a parallel implementation. We discuss the effect on the numerical stability and conclude that the consequences are only of minor importance as lon...
متن کاملCalculs pour les matrices denses : coût de communication et stabilité numérique. (Dense matrix computations : communication cost and numerical stability)
This dissertation focuses on a widely used linear algebra kernel to solve linear systems, that is the LU decomposition. Usually, to perform such a computation one uses the Gaussian elimination with partial pivoting (GEPP). The backward stability of GEPP depends on a quantity which is referred to as the growth factor, it is known that in general GEPP leads to modest element growth in practice. H...
متن کاملOn the Parallel Complexity of Gaussian Elimination with Pivoting
Consider the Gaussian Elimination algorithm with the well-known Partial Pivoting strategy for improving numerical stability (GEPP). Vavasis proved that the problem of determining the pivot sequence used by GEPP is log space-complete for P, and thus inherently sequential. Assuming P 6 = NC, we prove here that either the latter problem cannot be solved in parallel time O(N 1=2?) or all the proble...
متن کاملScalable Stable Solvers for Non-symmetric Narrow-Banded Linear Systems
Banded linear systems with large bandwidths can be solved by similar methods as full linear systems. In particular, parallel algorithms based on torus-wrap mapping and Gaussian elimination with partial pivoting have been used with success. These algorithms are not suitable, however , if the bandwidth is small, say, between 1 and 100. As the bandwidth limits the amount of parallelism available a...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2002